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Abstract. 



5— i ' We consider discrete spectra of bound states for non-relativistic motion in attractive 

potentials V a {x) = — \V \ |x| _cr , < a < 2. For these potentials the quasiclassical 
approximation for n — > oo predicts quantized energy levels e a (n) of a bounded 
spectrum varying as e a (n) ^ — n - 2a '/( 2 -< J )_ w e construct collective quantum states 
using the set of wavcfunctions of the discrete spectrum assuming this asymptotic 
behaviour. We give examples of states that are normalizable and satisfy the resolution 
of unity, using explicit positive functions. These arc coherent states in the sense of 
Klauder and their completeness is achieved via exact solutions of Hausdorff moment 
problems, obtained by combining Laplace and Mellin transform methods. For a in 
the range < a < 2/3 we present exact implementations of such states for the 
parametrization a = 2(k — l)/(3k — 1), with k and I positive integers satisfying k > I. 
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1. Introduction 

The construction of collective quantum states characterizing the whole spectrum of a 
quantum system is a challenging problem which depends sensitively on the nature of 
the potential involved. The standard coherent states (CS) are closely related to the 
harmonic oscillator potential and are defined, for complex z, by 

oo n 

|,)=exp(-|,| 2 /2)^^=|n), (1) 

n V m 

where {|^)}^Lo * s the F°ck space of eigenfunctions of operator H = p 2 /(2m) + mu 2 x 2 /2 
satisfying H\n) = hioin + l/2)|n), (n\n') = 5 njn >. There have been many attempts to 
generalize the construction of Eq. ([T]) for potentials other than ~ x 2 [TJ [2] . These efforts 
were hampered by the fact that the exact eigenstates and spectra of general potentials 
V(x) are known in only a few cases, and for a very special form of V(x). In fact for purely 
power-law type potentials of the form V a (x) = ±|Vo| \x\~ a the only exact solutions we 
know are for \o~\ = 1,2. The attractive case a > 2 is considered in a certain sense as 
unphysical j3] . This paucity forces one either to resort to approximations or to construct 
trial generalizations of (pQ) in which some known ingredients are built in, whereas other 
features are not accounted for. In this context, the quasiclassical approach may play a 
prominent role [U El [7] . 

In this work we will be concerned with attractive power-law potentials and the 
discrete part of their spectrum, which is known to be bounded [3] and which we assume 
here to be nondegenerate. Here V a {x) = —\Vo\ \%\~ a , < a < 2 and the quasiclassical 
estimate for the spectrum E a {n) is obtained from the Bohr- Sommerf eld quantization 
rule 

l - J v^m [E a {n) - V a {x)} dx = n(n + 1/2), n = 0, 1, . . . (2) 

which can be evaluated as (with H = 1) 

/ I 1 In \ -2a/(2-cr) 

EM = - ( \ -^i^T I Vol ) ^ -n"-/^), n -► oo, (3) 

see [6] for derivation and |11 15] for various refinements. In Eq. (j2j) x = a and x = b are 
turning points of the potential, defined by p(a) = p{b) = 0, and D(o~) = f a y/x a — ldx. 
We now use dimensionless units in which H a (p,x) =>- h a (p,x) = p 2 — x~ a ' . Then, 
we fix the value e a (0) = and thus arrive at the quasiclassical form of the spectrum 

e a {n) « l-cn~ 2a/( - 2 ~ a) + ..., < a < 2, n -> oo, (4) 

where the constant c > 0. We shall now follow the approach developed in [H HI [10] 
and incorporate the form of Eq. (J3J) to construct the generalization of states defined 
in Eq. ([1]) which are specially adapted to attractive potentials proportional to — |x| _CT . 
(For a recent treatment of CS for continuous spectra, see [HI ITT].) 
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To do so we use the method elaborated for discrete spectra and propose the specific 
form of a collective quantum state spanned by a set of eigenfunctions \n, a) of h a (p, x) 
satisfying (with (n, a\n',a) = 8 n>n >) asymptotically, as n — > oo 

h a (p, x)\n, a) = e a (n)\n,a). (5) 

Note that the equality sign in Eq. ([5]) is in general not valid for small values of n. 
Although \n,a) are not known in general, we still construct a trial wave function 
[S H2l [13] in the form generalizing Eq. ([I)) : 

I J, 7, cr) = N-nJ) £ J " /2eXP ^ (r2)) In, a), (6) 

where J > 0, and 7 are real and p CT (n) are so chosen as to assure the convergence of the 
normalization A r cr ( J) > 0, 

00 jn 

M a {J) = V —— < 00, < J < R < ex). (7) 

n=0 

The choice of p a (ri) in Eq. ([6]) is dictated by an "action identity" of the form [8], [9] 

(J, 7, a\h a {p,x)\J, 7,0-) = J (8) 

which directly implies 

n 
pAn) = Y[e a (j), p ff (0) = l, (9) 

and consequently the basic relation follows: 

e ff (n)= (\ v n = l,2,..., (10) 

p ff (n- 1) 

e ff (0) = 0, 

which, within our approach should be understood in the asymptotic sense. The states 
I J, 7,0") should satisfy the resolution of identity with a weight function W a (J) > 0: 

/oo 
dJd7 |J,7,o-)W (T (J){J,7,a-| = J^|n,<r)<n,o-| = I (11) 

n=0 

which reduces (wde Eq. (102) of Ref. [9]) to an infinite set of integral equations for an 
unknown positive function W a (x): 



R 

n 



x 





WJx) 



_N a {x)_ 



r R — 

dx= x n W a (x)dx = pa(n), n = 0,l,... (12) 

Jo 

where p a (0) = 1. If R < 00, Eqs. (Tl2|) is the Hausdorff moment problem [13]. It is 

known that if for a given set of p(n)'s the positive solution of Eqs. ( 112 j) exists then it 

is always unique [13] • The situation is very different for Hamiltonians with unbounded 

discrete spectra which lead to the Stieltjes moment problem with R = 00 in Eq. ( TP21 . 

In this case the solutions can be either unique or non-unique, see [15]. Observe that 



Quasiclassical Asymptotics and Coherent States for Bounded Discrete Spectra 4 

apart from their orthogonality no specific knowledge of the \n, cr)'s is required to derive 
Eq. (H2D. 

The state \J, 7, a) defined by Eq. (|6]), with p a (n) satisfying Eqs. ( ITOj) and ( 1T2l) . 
is a generalized CS state in a sense of Klauder [9], asymptotically relevant for 
V a {x) rsj - \ x \- a . 

2. Generating solutions of Hausdorff moment problems 

Our strategy from now on is: a) identify the form of p cr (n)'s to assure that Eq. ( 1T2|) can 
be solved for W a (x) > 0; b) calculate the associated energy spectrum from Eq. ( TTUj) ; 
and c) identify the exponent a obtained from Eq. (TJ2) and thus link the potential V a {x) 
to p a {n) and |J, 7, a). Evidently the correspondence in c) above is not unique: one may 
give different pi (n), r = 1, 2, . . . which yield the same asymptotics via Eq. (j4]). 

We consider Eq. ( IT2"1) as a Mellin transform /*(s) = Ai[f(x); s] = f Q °° x s ~ 1 f(x)dx, 
s complex [16]: 



W CT (x);s 



p a (s-l), s>l, (13) 



M 

or equivalently 

W a (x)=M- 1 [p a (s-l);x] J (14) 

where A^" 1 is the inverse Mellin transform. We observe that the moments used in this 
work will always be of such a nature that the integration range R in Eqs. (!T2|) will be 
equal to 1. Then the moment sequences will be decreasing functions of n. (We stress 
that this is not a general rule: there exist Hausdorff moment problems necessitating 
R > 1 for which the moment sequences are increasing [T7l HE].) 

In our search for positive solutions of Eqs. ( TT2|) we were greatly helped by a relation 
between the Laplace transform and a special case of Mellin transform [20]. To elucidate 
this link suppose that a function F(x) is considered for which its Laplace transform is 
known: 



C[F(x);p] = F(p)= e~ px F(x) dx, p > 0. (15) 

Jo 

We now perform a change of variable x = \n(l/y) in Eq. (fl~5l) . which gives 

POO 

F(p) = / y^ 1 F(\n(l/y)) H{1 - y) dy, (16) 



JO 

where H(z) is the Heaviside function. Through a formal renaming p -h- s we treat 
Eq. (fTBl as a Mellin transform 

F(s)=Af[F(ln(l/x))#(l-x);s]. (17) 

The relations of Eqs. (fl~5l)-( TT7|) allow one to search for possible solutions of the Hausdorff 
moment problem ( TT2l) for R = 1 via the method of the inverse Laplace transform 
[T9] . with a succession of following steps: i) choose a strictly decreasing sequence 
of moments p(n), n = 0,1,...; ii) rename them as F(p + 1), s.Eq. ( TTBT) : iii) search 
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for the inverse Laplace transform F(x) corresponding to F(p) and check whether 
W(x) = F(\ja(l/x)) H (1 — x) is a positive function on [0,1]; then, if so, W(x) is the 
solution of the Hausdorff moment problem Eqs. ( TT2l) . One is helped here by the fact if 
that F(x) is positive on [0, oo) then F(ln(l/x)) is positive on [0, 1]. 

2.1. Illustrative example 

We illustrate this approach with an example directly related to our construction. We 
choose the moments as p{n) = eexp(— \Jn + 1), p(0) = 1; the relabelling p{n) -B- F(p+1) 
gives F(p) = eexp(— y/p) which, with the formula 2.2.1.9 on p. 52 of [19] for a = 1, is 
the Laplace transform 



poo 

eexp(-y/p) = / e" 
Jo 



-px 



exp(— l/4x) 



dx. (11 



l\/ : Kx' i / 2 

In the next step we verify that j^[ln(l / x)]~ 3 / 2 exp I — j^kj^. ) is a positive function on 
[0, 1] and consequently 



x n 



e 

exp 



L2 v ^F[ln(l/a:)]3/2 ~* \ Un(l/x)J\ ** ~ 6 ' 6 ~ " ' (19) 

n = 0, 1, . . ., is a complete solution of the Hausdorff moment problem Eqs. (1T21) . With 
the above moments the spectrum Eq. ffTUl) is 

e( n ) = e^-v^+T, (20) 

and its n — > oo asymptotics is 

e(n)^l--\- + ^ + ..., (21) 

which with Eq. (j4j) determines o = 2/5 and c = 1/2. The generalized coherent state 
describing such a spectrum is given by Eq. ((6]) with the normalization 

1 °° 
A/"( J) = - J^ e (n+1)1/2 J n , < J < 1. (22) 

n=0 

The CS | J, 7,2/5) is then asymptotically relevant for motion in the potential V2/s(x) ~ 
_| a .|-2/5_ fj^g f orrnu i a Eq_ f fT9|) can also be cross-checked by referring to the tables of 
inverse Mellin transforms (see formula 3.7 on p. 174 for a = 1 of Ref. [20]). In the 
spirit of Eqs. ( fi~2l) we call the weight function in Eq. ( |T9l M^2/s(x), < x < 1, which 
can also be derived from eL(l/2,ln(l/x)), where £(7, a;) = \^ x ~ 3 ^ 2 e~" / ^ 2x " > is the 
so-called one-sided Levy stable distribution [2T] . 

We now present a more general case which can be treated with the above method. 
For that purpose we again stress that we are looking for special sequences of moments 
p{n) satisfying the Hausdorff moment equations which at the same time possess very 
specific asymptotic properties implied by Eqs. (@J and ffTUj) . These are very restrictive 
conditions indeed. The search for such solutions may proceed by exclusion and at the 
beginning it was not certain at all if a general solution existed. It is all the more satisfying 
that a parametrization can be given that produces a vast ensemble of solutions, at least 
for some range of values of a. 
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2.2. Full solutions for a rational in the range < a < 2/3 and for a = 1 

Let us define a sequence of moments, parametrized by a, u, k and I, and given by: 



p M (k,l;n) = e a {n + l)- v e- a{n+l)/ , n = 0,l,... 



(23) 



with conditions: k and I positive integers; k > I; a > and v > 0. The last condition 
assures the positivity of the weight function, see Appendix A. We use now the formula 
2.2.1.19 listed without proof on p. 53 of \19\: 

Vkl 1/2 ~ u „_i^ n fa k l l 



-px 



(2tt)(*-0/2 






A(fc,0) 



cix = p-^ e - a 



,!/fc 



(24) 



for p > 0, where G™' n (z\ • ■ ■) is Meijer's G function [22]. The detailed demonstration of 
Eq. f[2~4"]) will be given elsewhere. 

We transform the Eq. (124)) with x = \n(l/y) and arrive at the expression of the 
type of Eq. ( TIB]) , namely 

(25) 



.P- 1 



(27r)(*-0/2 



x g; 



/,;.() 



fc ;Z 



a fe Z 



Z 1 / 2 - [^(l/y)]- 1 



aK) H{1 - y) 



dy 



e a p- u e- ap 



l/k 



p>0. 



< fc \k k [\n(l/y)} 1 
In Eqs. (J24l) and (l25l) we use a compact notation for special lists of k elements |19j : 



A(M) = f,^, 
transform 



P>9 



a+ k 1 . The Meijer G function is defined as an inverse Mellin 



ai...a p 
01- A 



M 



-i 



n?= m+ ir(i-/3 J - s )n^„ + i r K + s ) 



(26) 



= G([[«i, . . . , a n ], [a n+ i, . . . , a p ]], [[A, . . . , /3 m ], [/3 m+ i, . . . , /3 g ]], z), (27) 

where in Eq. (126]) empty products are taken to be equal to one. In Eqs. (1261) and (J2j 
the parameters are subject of conditions: 

z 7^ 0, < m < g, < n < p; 



(2? 



«j£ 



j = l,...,p; /3j G C, j = l,...,q. 

For a full description of integration contours in Eq. (126]) . general properties and special 
cases of the G functions see [25] [23]. In Eq. (I2T]) we present a transparent notation, 
which we will use henceforth, inspired by computer algebra [23]. With this notation 
Eq. (I25p now becomes 

nmr 1 (29) 



y 



p-i 



(27r)(*-0/2 



x G [[], [A(Z, ,/)]], [[A(fc,0)],[]], 



fc 7« 



a fe Z 



fc fc [ln(l/y)]' 



dy 
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= e a p- u e~ apl/ \ p>0, 



or, when rewritten as the moment problem with Eq. (T53 



i 

y n 





e a Vkl 1 / 2 ~ u ri .„ . ..„ , „ /.. . rA /7 ... rrA „ „,. r nn a fc / 



k H 



(27r)(fc -0/ 2 HlM^G^l ],[A(/,,)]],[[A(fc,0)],[]] i fcfcNlM J 
p M {k,l;n) (30) 

^M^(M;</)^, « = 0,l,...,oo. (31) 

Observe that in Eq. (j30|) only the second and third lists of parameters are non 
empty in the G functions (vide the notation of Eq. (1271) ). as may be inferred from the 
conditions of Eq. ( |28l) . Eqs. ( 1301) and (13~TT) have a very rich ensemble of solutions which 
will be studied for various values of the parameters a, u, k and I. The role played by a 
and v is fundamentally different from that played by k and /. The asymptotic behaviour 
as n — > oo of p^ a ^ (k, I; n) does not depend on a and v and the exponent of the power-law 
n dependence of spectra is a function of l/k only: 

e ^) (fcj/;n) = i_^_^_ + ... j (32) 

which, by Eqs. ([3]) and (J3J) immediately implies that a may be "fine-tuned" with the 

parametrization 

2(k — I) 
*= 3 fc _/ , fc>^5 A;,Z = 1,2, ... . (33) 

Eq. fl33|) confines a to a possible range of < a < 2/3. In Eq. ( 132]) the constant c of 
Eq. 01]) is equal to c = al/k. The above analysis indicates that the weights W^ a ' v ^(k, I; y) 
for different a and v give the same asymptotics. 

We shall give examples of such an asymptotic " degeneracy" with explicit forms of 
W^ a ' v \k, I; y) 7^ W^ a ' u \k, I; y) for a few v ^ iS '. For simplicity, from now on a fixed 
value a = 1 will be used for all examples derived from Eq. (I24J1 . We shall observe the 
onset of complexity with increasing values of k and /: starting with k = 4 and I = 1 
we leave the realm of standard special functions as then the corresponding Meijer's G 
functions can be only converted to finite sums of generalized hypergeometric functions 
of type pF q . They are however available through computer algebra systems [23] and 
their properties are, to a large extend, readily accessible. Since a = 1 we shall denote 
W M (k,l;y) = WM(k,l;y). 

2.3. Special cases: 

1. k = 2, / = 1 and u = 0: 

eU/ii 



W (0) (2,i;y) = 3rT^G [[],[o]], [[0,1/2], [ ]], 



ln(l/y) V L J ' L JJ ' U ' ' J ' LJJ '41n(l/j/) 



G f [[],[]], [[1/2], []], -n4-7-v )■ 



ln(l/y) V U m JJ ' U ' J ' LJJ '41n(l/y)y 
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which, as expected, reproduces precisely Eq. (fT9"|) . 
2. k = 3, I = 1 and u = 0: 

^ (0) (3, 1; 2/) = ^^ G ([[ ], [0]], [[0, 1/3, 2/3], [ ]], ^^ ) (35) 
= wK [[U]U[1/3 ' 2/3U]] '27l^) 



— [\n(l/y)r 3/2 K l/3 f^M., 
3tt [ V /y;j 1/3 I 9[ln(l/?/)] 1 / 



[ln(l/y)]^ 3 / 2 Ji 1/3 ( nnTiM V /a ,, 11/2 ) , (36) 



where K v {z) is the modified Bessel function of the second kind. Eq. (1361) can also be 
checked with formula 3.13, p. 175 of [20]. With Eq. (132]) the corresponding spectrum 
varies as 

e( 1 '°>(3,l;n)-H-^ i7 3 + ..., twoc, (37) 

yielding a = 1/2 and c = 1/3. 

3. k = 3, Z = 1 and v = 1/2: 

^ (I/2> ( 3 '^» = R^ G (tn.[^])J[0,l/3,2/3]J]],^ IM 

= 5(^37? I'"' 1 /!/)]- 1 *V» ( 9|lll( ^ )|1/2 ) ' *V3 ( 9[ln( ^ )]1/2 j • (38) 

which yields the same n- dependence as in Eq. (|37|) . 

It is instructive to compare weight functions leading to the same spectrum 
asymptotics. To this end we present the weight functions from Eqs. fl37|) and (138|) 
in Fig. Q] . 

4. fc = 3, Z = 2 and i/ = 0: 



^ (0) ( 3 '2;y) = ^77^G( ], [1/2]], [[1/3, 2/3], [ 



ln(l/y) V UJ ' L ' JJ ' LL/ ' ' J ' LJJ '27[ln(l/y)P 
which can be neatly expressed in terms of modified Bessel functions K v (z): 

^>(3, 2 ;,)^HiM- 3 ex P (^-i 7 - F ) 



X 



% L7,J,uJ +%' 



27pn(l/y)]V z/ ° V27[ln(l/y)] 



12 



(39) 



which leads to the asymptotics: 



e (^)(3,2;n)^l-3^ + ..., (40) 



with a = 2/7 and c = 2/3. 
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W(y) 




Figure 1. Plot of the weight functions W(y); the line / corresponds to W^°'(3, l;y), 
see Eq. (|56| and the line II represents W^^ 2 \3, 1; y), see Eq. 



5. k = 3, I = 2, v = 1/4: 

This case can be expressed by the G function 

2- 1 /*V3 



W^ 4 \3,2;y) 



[ln(l/y)]3/4 



xG {[[], [1/8, 5/8]], [[0,1/3, 2/3], [ ]], 



J '27[ln(l/y)]V' (41) 

which has a representation in terms of a sum of three hypergeometric functions of 
type 2-^2) which will not be quoted here. The asymptotics is that of Eq. (T4"0"j) . The 
weight functions from Eqs. ( 139]) and (j4l~j) which share the same spectrum asymptotics 
are compared in Fig. |2j 
6. k = 4,l = l,u = Q: 

The corresponding weight function W^(A, 1; y) has an exact representation in terms of 
a sum of three hypergeometric functions of type 0-P2 which we will not given here. 

This pattern extends to higher values of k and I for which the weight functions 
become increasingly complicated. They can however be fully handled analytically and 
graphically with a relative ease: we always deal with finite sums of hypergeometric 
functions. 

We shall go over to further examples and employ a wealth of formulae of Ref. [19] 
and [20] , different from that of Eq. (E 



7. We now choose a function differing from the exponential which nevertheless yields 
asymptotic behaviour which is close to Eq. (l2il : this is provided by the formula 3.16.6.7, 
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W(y) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

y 



Figure 2. Plot of the weight functions W(y); the line / corresponds to W^°'(3,2; y), 
see Eq. <J3H]) and the line II represents VK (1/4) (3, 2; y), see Eq. (|2Tj) . 



p. 358 of [19] or, alternatively by Eq. 7.69, p. 230 (for a = 1) of [20], namely: 



2#„(1) v/7rln(l/y) 



exp(- 



■ ln(l/j/)- 



AT, 



i//2 






iln(i/y)>/ 



r/,(/ 



i ^(v^+T; 



/#<>), n = 0,l,.... 



(42) 



K V {1) v^TI 

The modified Bessel functions K u (z) in Eq. (142]) for z/ 7^ p/2, p = 1,2,... are not 
elementary functions. The weight function in the l.h.s. of Eq. f !42p is positive on [0, 1] 
and normalized as pi (0) = 1. The asymptotic behaviour for v = 4/3 is close to that 
of Eq. (J2U): 

1 2 -, (43) 



e (K) (n) ~ 



2nV2 8n 
leading to c = 1/2 and o = 2/5. 

The weight functions from Eqs. (fl~9l or (|35l) . and Eq. (144]) display the same spectrum 
asymptotics and are illustrated in Fig. |3] 

8. We use now Eq. 2.2.2.1, p. 53 of [19] for the choice v = 4/3 and a = 1, which leads 
to the normalized Hausdorff moment problem: 



e- 1 [ln(l/y)] 1/6 / 1/3 (2VlnTl7^) 



dy 



,l/(n+l) 



in 



!)4/3= 



giving the asymptotics with a = 2/3 and c = 4/3: 

4 
e(n) ~ 1 - — + ... . 
on 

In Eq. ( 145]) I v {z) is the modified Bessel function of the first kind. Note that the value 
a = 2/3 cannot be obtained from Eq. 



(44) 



(45) 
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W(y) 




Figure 3. Plot of the weight functions W(y): the line / corresponds to W^°'(2, l;j/), 
sec Eqs. (|3~i)) and (|TTJ)) and the line II represents the weight function in Eq. (j^)) . 



9. In this final example we address the problem of the one-dimensional Coulomb 
potential ~ — |x| _1 for which the exact spectrum is [3] 

1 



e c {n) 



[n 



\2- 



n 



0,1, 



oo. 



(46) 



The CS for this case have been constructed in [8] and the corresponding moments 



n+2 



n 



are given by p c (n) - - 2(n+1) , 
corresponding weight function 



0,1, 



and they lead to the exact form for the 



I")] H(l-y). 



(47) 



The Coulomb problem, even in its simplified version treated here, presents a particularity 
in that its resolution of unity is distributional in character as it involves, in Eq. ( 1471) . 
the Dirac delta function S(z). Although the solution Eq. (|47|) is unique for the exact 
moments p c {n) defined above, we are able to construct another weight function V c (y) 
which will asymptotically reproduce e c (n) of Eq. ( 1461) . To this end we use the formula 
2.2.2.8, p. 54 of pi] for a = 1, which gives 



V n V c (y) dy 



V" 



h{2yjMTJv)) 

[ln(l/y)]V2 



S(y 



dy 



e -l e l/(n+D. 



n 



0,1, 



(41 



with asymptotics e c (n) —$■ 1 — n~ 2 + 3n~ 3 + . . . . Note that V c (y) displays a non-trivial 
dependence on y G [0, 1] but still retains a Dirac peak at y — 1. We are not aware of 
existence of any weight function without Dirac's delta leading to Eq. (14"6"|) . 

The different weight functions relevant for the Coulomb interaction, both of them 
involving Dirac's delta function at y = 1, are schematically displayed on Fig. SJ 
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W(y) 



0.2 



0.4 



0.6 



Figure 4. Plot of the weight functions W(y): the line I corresponds to W c (y), see 
Eqs. (|47p and the line II represents the weight function V c (y), see Eq. f4"5|) . The delta 
peak at y = 1 are represented by a vertical line. Observe that this is a schematic plot: 
different coefficients (1/2 and e _1 . s. Eq. (|4"T|) and P5|) respectively) in front of delta 
peaks at y = 1 cannot be accounted for. 



3. Discussion and Conclusion 



In this work we used a semiclassical approximation for the spectra of one- dimensional 
inverse power-law potentials to construct approximate coherent states, relevant for these 
potentials. Our goal in using this construction was to achieve the resolution of unity 
for a given potential characterized by an exponent a in V a (x) ~ — \x\~ a . In this sense 
the po-(n)'s of Eq. (Q should be perceived as a sort of trial parameters which should 
asymptoticaly reproduce e a {n) via Eq. ( ITUl) . This leads, of course, to a multiplicity of 
possible choices of p a (n) leading to the same e a (n). The price of this approximation 
is the fact, that the temporal stability characterizing exact Gazeau-Klauder CS [8j [TTj 
denoted by \J, 7, cf)gk (i-e. the equality e~ lh,Tt \ J, 7, <j)gk = \J,1 + t,°~)GK) will be only 
approximately satisfied by the states of Eq. (JS]). On the contrary, the basic Gazeau- 
Klauder axiom of the resolution of unity is fully maintained in our construction. 

For certain values of a, more precisely for rational a such that < a < 2/3 and for 
a = 1, we are able to produce many resolutions of unity that are asymptotically relevant 
for one and the same a. This is due in the first place to two key formulae, Eq. (j24|) 
and ( |33l) which involve the use of the inverse Laplace transform. For other values of a 
in the range 2/3 < a < 2 (except a = 1) our approach does not produce any required 
solutions for the resolution of unity. 

We should comment here on the nature of approximation involved in the formulation 
of Eq. ([6]). Since for general a and n neither exact spectra nor exact eigenf unctions 
I n, a) are known we strike a compromise in Eq. ([6]) by retaining the exact orthonormal 
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eigenfunctions and replacing the energy spectrum by its quasi-classical form of Eq. 01]). 
In Eq. (jSJ) we substitute the quasiclassical approximation for e a (n) instead of the 
exact spectrum. This means that we are using the asymptotic approximation in the low 
energy region, where its use is not a priori justified. However, for certain values of a 
the quasiclassical approximation is very successful even down to the low energy: in fact, 
for the linear repulsive potential V(x) ~ x [x > 0; V(x) = oo, x < 0) it predicts the 
correct energies and wave functions right down to the ground state [HI IE]- For attractive 
potentials, for which in general one does not have exact solutions, the agreement is less 
spectacular [6] but in general the quasiclassical approximation works well in large parts 
of the spectra and not only for n — > oo, in which region it tends to the exact solution. 
Therefore we believe that the use of Eq. (J4J) in Eq. ([6]) is a reasonable prescription. 
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Appendix A. 

We give here a streamlined proof of positivity W^' u \k, I; y), for < y < 1 from Eq. (1301) 
under the conditions k, I positive integers, k > I, and v > 0. 
Consider first Meijer's G function of Eq. ( 124"|) . with A(k,a) 



a a+1 

fc' k ' 



q+fc-1. 



G([[],[A(l,»)]],[[A(k,0)],l]],z)=M- 



n 



£r(« + £ 



r l-l 



u+j- 



n;=or(* + ^ 



-\z 



M' 1 



jj r(* + | 



fc-i 



nr< s+ 



J-J- T(s + £±21 
Vi=0 L ^ + I > / 3=1 

The convolution property for JH[f(x); s] = f*(s) and Ai[g(x); s] = g*(s) 

r°° dt r°° dt 

M- 1 [r(s)g*(s);x}= f( x / t )g(t)-= g(x/t)f(t)- 

Jo l Jo t 

for f(x) > and g(x) > clearly conserves the positivity, see [15] and references therein. 

Fix v > and consider A4 _1 of an individual term of the first product in (IA.2J) . It will 

obey, using the formula 8.4.2.3, p. 631 of [23] 



(A.l) 
(A.2) 

(A.3) 



AT 1 



T(s + ^ 



x 



x- 



■j/k M 



X 



,-i+iki¥l 



j=0,l,...,/-l,(A.4) 



r (H- + i¥)) 

which is a positive function for < x < 1. Also, concerning the second product in (1A.2I) 
an individual term evidently gives 



M 



r|s + t 



x 



x j/k e 



> 0, x > 0, j = /,..., fc-1. (A.5) 
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Then ( 1A.2J) can be viewed as a multiple, fc-fold Mellin convolution of positive functions, 
which by (1A.3J) is itself positive. 

The proof of positivity of W^ 1 '^ (y) is completed by remarking that if F(y) is positive 

on [0, oo) then for a and j3 arbitrary In ( - ) F I In ( - ) J is positive on [0, 1]. We 

stress that the v < case destroys the positivity of W^ 1,v \y) and is not relevant for our 
purposes. 
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